Ultra-High Energy Quenching of the LPM Effect: Imphcations for GZK- Violating 

Events 



(N 
O 

o 



> ' 

in 
in 

o : 

(N 

O ' 

■ 

6 



John P. Ralston^, Soebur Razzaque^ and Pankaj Jain^ 
^Department of Physics and Astronomy, University of Kansas, Lawrence, Kansas 66045 
^Department of Astronomy & Astrophysics, Penn State University, University Park, Pennsylvania 16802 ana 

^Physics Department, I.I.T. Kanpur, India 208016 
(Dated: February 2, 2008) 

Cosmic rays in the regime of lO^'' — 10^" eV should not exist, according to the GZK bound. 
Conflicts between data and the bound have led to suggestions that laws of physics must be modified. 
We find that physical processes neglected in the early history terminate the LPM effect, important 
in calibrating the energy of electromagnetic and hadronic showers in the GZK energy regime. The 
processes neglected are the direct losses of electrons and positrons, overlooked by LPM and greatly 
augmented by ultra-high energy circumstances. We present an exact formula for direct energy losses 
of e'^e~ pairs in electromagnetic showers. Numerical integration in current hadronic models shows 
that actual energy losses are always substantially larger than LPM estimates, but the energy at 
which crossover occurs depends on models. 

PACS numbers: 12.20.Ds, 14.60.Cd, 14.70.Bh, 41.60.-m, 95.30.Cq, 98.70.Sa 



Some thirty years ago, Greissen, Zatsepin and Kuzmin 
(GZK) [| found that cosmic rays above about 5 x 10^^ cV 
should not be observed. Yet since its inception the GZK 
bound has never been confirmed, as evidenced by obser- 
vation of dozens, if not hundreds of anomalous events. 
Extraordinary contradictions are indicated by angular 
correlations |g[ of GZK-violating showers and cosmolog- 
ically distant sources. Experimental facilities such as 
HIRES and AUGER are poised to explore the question 
with vastly more data. The conflict between the data 
and the bound has led many theoretical workers to sug- 
gest that laws of physics must be modified. Sugges- 
tions of new physics include exotic particles, magnetic 
monopoles [Q, neutrinos with large cross sections from 
low-scale gravity |^], and even the rejection of Lorentz 
invariance Q. 

High energy air showers tend to be measured better 
with increasing energy, due to increased signal-to-noise. 
However the calibration of shower energies depends on 
details of the shower evolution. Hadronic showers have 
some model dependence, but the electromagnetic com- 
ponents are thought to be well-understood. Many years 
ago. Landau and Pomeranchuck (LP) Q showed that 
high energy electromagnetic interactions in matter de- 
structively interfere, effectively increasing the interac- 
tion length. Migdal (M) gave formulas smoothly in- 
terpolating between different coherence regimes. Monte 
Carlo codes incorporate these effects, currently believed 
to extend to asymptotically high energies. The question 
of ultra-high energy calibration involves what and how 
much to believe from LPM. 

The LPM calculations hinge on dominance by a con- 
cept called the "formation length" of bremsstrahlung. 
In a peculiar remark, Galitsky and Guerevich (GG) ||] 
noted an instance where the absorption length might be 
less than the coherence length. Then photons would be 



absorbed before they were fully emitted. GG were dis- 
cussing e+e" pair production via the Bcthe-Heitler (BH) 
cross section. Apparently GG missed the suppression of 
pair absorption by LPM effects, which decreases the pair- 
production rate self-consistently. In fact the LPM regime 
extends to arbitrarily high energies, if treated as an elec- 
trodynamics problem in isolation. Yet nobody working 
before 1960 could anticipate the evolution of high energy 
hadronic physics. 

We highlight a rather subtle point: the hadronic com- 
ponent of electrons, positrons and photons. High en- 
ergy electrons undergo two independent causes of electro- 
magnetic energy loss: the radiative losses and the direct 
losses. Radiative losses generate amplitudes proportional 
to the acceleration (3. Direct losses exist as /3 0, and 
were not mentioned by LPM. Direct losses traditionally 
start with the equivalent photon density dN^/dx-y per 
Feynman Xj in the uniformly moving electron's fields. 
When such photons interact with cross section (Jtot , then 
the equivalent energy loss per unit distance z is given by 
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The electron then has hadronic interactions because atot 
has hadronic contributions, commonly called photonu- 
clear reactions. That is, the physical electron cannot 
be separated from a comoving cloud of photons, and 
the physical comoving photons cannot be separated from 
their comoving component of hadrons. 

Numerous hadronic interactions are carefully imple- 
mented in current air-shower Monte carlo codes. For 
muons and tau leptons the direct photonuclear losses 
are larger than the radiative losses and are included 
via equivalent photons in current codes: but not for 
the electron. Meanwhile the photonuclear reactions are 
not LPM-suppressed because they occur on very short 
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hadronic distances and time scales. Absorption of pho- 
tons also affects the LPM radiative losses. Since electron 
and positron pairs, plus photons, are the primary com- 
ponent of electromagnetic shower cascades, the photonu- 
clear losses by electrons may have an observable effect in 
shower evolution. 

Shower energy calibration also involves detailed fea- 
tures of the detector response, which is not our task. Our 
task here is the termination of the LPM regime above a 
certain energy. The exact electromagnetic energy loss 
comes from the Foynting vector. 
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Here Euj{b), B,^{b) are the electric and magnetic fields in 
Fourier to space evaluated at transverse spatial coordinate 
h. The trajectory of a fast charge has a series of small 
angle kinks distributed along a nearly straight path. Up 
to a phase, the radiation fields from a particular kink are 
given by 
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where n is a unit vector pointing to the observation point 
and A/3 is the change in velocity at the kink in units of 
c. Kinks separated by less than the formation length of 
radiation add their fields with destructive interference. In 
free space the formation length Lq ~ 7^/^ for a charge 
with Lorentz boost factor 7, while in the LFM regime 

The kink field [Eq. (^] transforms like a boosted 
dipole. Meanwhile the "kinematic" monopole fields of 
the boosted Coulomb charge are azimuthally symmet- 
ric about the trajectory. By orthogonality over the az- 
imuthal integral, there is zero net interference between 
the energy losses of the boosted Coulomb fields and the 
radiation fields: 



dE ^-^radiativc ^ ^-^dircct 

dz 



dz 



dz 



(4) 



We set aside the radiative losses for now, and turn to 
the direct term. The fields of a uniformly moving charge 
are given by 
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where Kq and Ki are modified Bessel functions. The 
effective Lorentz boost factor in the medium, denoted by 
7ni, is defined as 
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The corresponding formula for the energy loss is 
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attributed to Fermi pO| . Here 6min is the minimum im- 
pact parameter for which the loss subprocess proceeds. 

The Fermi formula is rather more sophisticated than 
the equivalent photon approximation. Replacement of 
1/e — 1 ^ (e — 1) ^ cr plus the limit a; ^ 0, are needed 
to find the approximation Eq. (|^). Equivalent pho- 
tons can be used to calculate various subprocesses da, 
but only at leading-log order. When subprocesses are 
pursued beyond leading order certain re-summations are 
needed, which depend on the process. But there are no 
restrictions on Poynting's theorem, even in the quantum 
domain. All the terms needed for energy loss are summed 
in the "Fermi formula" Eq. (^, which is deeper than it 
may appear at first sight. 

As a consequence, there are two distinct regimes: 
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We use the "absorption length" 1 / (a; Im[e^]) = Labs as 
a characteristic length, of order l/natot- The regime of 
Eq. (H) corresponds to free space, or any other cases 
of negligible medium effects. The Bessel functions then 
permit lu < 7717, which is the kinematic maximum. Then 
dE/dz ~ E for the familiar case of Labs being con- 
stant. In the LPM model 1/Labs ~ V^/B, (jJmax is 
the kinematic maximum at high enough energies, and 
dE/dz ^ \fE . In contrast, the other regime of Eq. 

is dominated by absorption in the medium, and 
'^max ~ m/-\/|et^ — 1|. The corresponding loss integral is 
constant if (say) Labs is constant. The crossover between 
the two regimes occurs when the absorption maximum 
falls below the kinematic maximum, namely 
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absorption dominates. 



(10) 



Mutiplying both sides of Eq. ( |lo| ) by wavelength 
A, the regime of absorptive dominance corresponds to 
7^A > Labs- The left-hand side is the free-space forma- 
tion length, which ceases to apply when it exceeds the 
absorption length. This is just the criterion for absorp- 
tion dominance cited by order of magnitude formulas by 
GG §. 

We continue by noting that absorption and propaga- 
tion are linked by unitarity and dispersion relations relat- 
ing the real and imaginary parts of propagators. Indeed 
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the renormahzed photon propagator in a medium is el- 
egantly coded into the permittivity appearing in 7^ 
and Eq. (|7|). By definition 

1 1 
^ = = , (11) 

e^Lj'^ — fc2 Cij2 — fc^ — E 

where E is the irreducible self-energy from summing 
bubble diagrams. (The "pole model " for of non- 
relativistic tradition is one example.) Either symbol, 
or S is equally valid to parametrize the propagator, with 
the identity relating them e^j — 1 = In this form 

the optical theorem is 



lm[e^] = 



ncrtot(w) 
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where n is the number density of scatterers with pho- 
ton total cross section (Ttot(^)- Given crtot, then the full 
dielectric constant follows by dispersion relations. 



Re[£{uj)] = 1 + ^75 J 
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where V is the principal value. 

We turn to the total cross section o'tot(w). In the LPM 
model the pair production cross section is 
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where Eg = 21.2 MeV is the scattering energy scale and 
Xq is the radiation length. The corresponding dielectric 
properties are 
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where i?LPM is the LPM scale. Note that o-lpm(ijj) de- 
creases with UJ, as necessary for the self-consistency of 
the LPM argument. 

Meanwhile photonuclear contributions increase with oj. 
At some high energy, the LPM-model absorption cross 
section can only be a small fraction of Ctot- This is 
because hadronic total cross sections typically scale like 
tJ ~ s°, where s — 2mpE is the center of mass energy- 
squared on a target of mass mp. Two explanations for 
this dependence are Regge/Pomeron theory, and the ob- 
served growth of small-x parton distributions q{x). Data 



from HERA shows that xq{x) 
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to a good approx- 



imation, yielding crtot ~ E^ '^. Growth with energy is 
simply due to the increased number of partons, and oc- 
curs in just the same way for the neutrino-nucleon cross 
section ||l^. The 1998-2001 photon cross section by Don- 
nachie and Landshoff (DL) ||l4 is a highly cited example: 
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We will use DL as an up to date baseline for cr^, although 
all extrapolations into the UHE regime should be used 
with caution. The corresponding dielectric properties are 
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where o^min ~ GeV, F{a, (3, 7, 5) is the hypergeometric 
function, Cj are the prefactors and dj are the powers 
of ijj in Eq. (p^). Cross sections are expected to grow 
until parton saturation sets in, after which asymptotic 
behavior such as the Froissart bound is expected. The DL 
total cross section reaches a modest geometric value (40 
mb) a± E ^ 10^^ eV, and extrapolations far above that 
value might be slowed to log'^{E) or similar asymptotic 
behavior. 

We are not aware of a model for crtot {^) which does not 
rise in the high energy regime. We compared the 1981 
fit of Bezrukov and Bugaev (BB) such that Froissart 
bound dependence sets in at about 50 GeV. The BB fit 
is 
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We use the BB model to compare model-dependence. 
We found analytically by dispersion relations in terms 
of poly- log functions. A more useful representation as 
asymptotic series is: 
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Potential effects for GZK are illustrated in Fig. [|, 
based on numerical integration of Eq. (Q) for air at 
sea level {Xo ^ 30420 cm; bmin — l/n^7r)- Conventional 
radiative losses (also shown) reveal the LPM downturn 
above eI;^'^^ - 2.3 x lO" MeV. The DL-based calcu- 
lation becomes non-negligible above about 10^^ eV and 
exceeds the radiative losses above 5 x 10^" eV. The BB- 
based losses are negligible through this energy range, but 
also exceed the radiative losses above 10^^ eV. Eventu- 
ally any crtot (t^^) model on the market will terminate the 
LPM regime. 

The onset of large direct losses of the DL-model (Fig. 
|l|) coincides alarmingly with the energy range where the 
GZK bound is confronted. As far as known, the actual 
cross sections could easily be larger than DL's, caus- 
ing the regime to occur even earlier. However we have 
not finished with the repercussions, because there re- 
main the radiative losses postponed earlier. For any 
given = oj/E fraction, the LPM coherence length 
^LPM ~ increases indefinitely with 7. Meanwhile 
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FIG. 1: Total energy loss [dE/dz) of electrons in air at sea 
level {Xo = 30420 cm) assuming standard LPM radiative 
losses. Losses with the DL model are indicated by the solid 
line, with the BB model (long-dashed) shown for comparison. 
Short-dashed line shows the LPM-suppressed radiative losses, 
which turn over at i^LPM ~ 2.3 x 10^^ MeV. The direct losses 
of both the DL- and BB- models terminate the LPM regime 
(radiative losses) above 5 x 10^*^ and 10^"^ eV, respectively. 

the photonuclear absorption length Labs ^ (72^) ° de- 
creases. When Xlpm > ^abs, the radiative losses are no 
longer comoving and superposing "on top of one another" 
near the speed of light. Instead, the radiation from each 
kink is absorbed before it reaches the next kink. Thus 
quenching of radiative LPM occurs just where the di- 
rect losses become noticable This is just the same 
physics as the experimentally observe d low energy 
"polarization effect" of Ter-Mikaelian ]1CI|, postponed to 
a regime of large enough e — 1 . We have nevertheless pre- 
pared Fig. 1^ using the standard LPM radiative losses to 
facilitate comparison. At an order of magnitude energy 
(say) above the crossover point, the direct losses domi- 
nate anyway, so that it does not matter that the LPM 
formulas are an overestimate. The crossover from LPM 
radiative to direct absorptive regime is inevitable. Unfor- 
tunately the uncertainties of ctot make the onset energy 
of absorptive dominance very uncertain. 

Figure has implications for the calibration of air 
shower energies. The maximum number of shower par- 
ticles iVmax increases with primary energy. The effects 
of LPM are to slow this increase: for photon-initiated 
showers at 3 x 10^'' eV, A^^ax 1-3 x 10" with LPM 
suppression, and TVinax ~ 2.2 x 10" with LPM turned 
off |15{. Adding direct losses is akin to turning LPM off 
(Fig. 111). However observational constraints on shower 
energies include not only A'max, but also the shape of the 
shower based on Greissen-type shower profiles. There is 
no direct way to go from a particular cross section to the 
entire shower profile, which incorporates every possible 
particle reaction. Indeed, the direct losses we study ap- 
ply not only to the electron, but to every single charged 
particle! 

We conclude that extrapolations of energy calibra- 



tions to 10 eV and above are incomplete and model- 
dependent. Large photonuclear cross sections violating 
no existing physical principles may have drastic effects. 
The effects on observations remain to be determined, and 
likely cannot be resolved without a concensus among ob- 
servers. With hundreds of papers citing a looming GZK 
crisis, we find a new burden of proof must be borne by 
those citing violation of the GZK bound. Observations 
in the regime of 10^° eV will likely be controversial for 
some time to come. 

Work supported in part under Department of Energy 
grant. We thank Graham Wilson and Doug McKay 
for comments. SR was partially supported by NSF 
AST0098416 grant. 



[1] K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G. Zat- 
sepin and V. Kuzmin, Pis'ma Zh. Eksp. Teor. Giz. 4, 
114 (1966) [JETP Lett. 4, 78 (1966)]. 

[2] D. Bird et al, Astrophys. J. 441, 144 (1995); J. El- 
bert and P. Sommers, Astrophys. J. 441, 151 (1995); 
N. Hayashida et al, Astrophys. J. 522, 225 (1999). 

[3] G. Farrar and P. Biermann, Phys. Rev. Lett. 81, 3579 
(1998); Y. Uchihori et al, Astropart. Phys. 13, 151 
(2000).; A. Virmani et al, Astropart. Phys. 17, 489 
(2002); G. Sigl et al, Phys. Rev. D 63, 081302 (2001). 

[4] T. Kephart and T. Weiler, Astropart. Phys. 4, 271 (1996) 

[5] P. Jain et al, Phys. Lett. B 484, 267 (2000); A. Jain et 
al, Int. J. Mod. Phys. A 17, 533 (2002). 

[6] L. Gonzalez-Mestres, in Proc. of the 25th ICRC, 1997, 
Durban, ed. by M. Potgieter et al (World Scientific, Sin- 
gapore), 6, 113; G. Amelino-Camelia et al, Nature, 393, 
763 (1998); R. Gambini and J. Pullin, Phys. Rev. D 59, 
124021 (1999); S. Coleman and S. Glashow, Phys. Rev. D 
59, 116008 (1999); F. Stecker and S. Gla show, Astropart 
Phys. 16, 97 (2001); J. Carmona et al, ^iep-th/020715^ . 

[7] L. Landau and I. Pomeranchuk, Dokl. Akad. Nauk SSSR 
92, 535 (1953); ihid 92, 735 (1953); A. Migdal, Phys. 
Rev. 103, 1811 (1956); R. Blankenbecler and S. DreU, 
Phys. Rev. D 53, 6265 (1996). 

[8] V. Galitsky and I. Gurevich, Nuov. Cim. 32, 396 (1964). 

[9] S. Dutta et al, Phys. Rev. D63, 094020 (2001). 
[10] E. Fermi, Phys. Rev. 57, 485 (1940); M. Ter-Mikaehan, 
High-Energy Electromagnetic Processes in Condensed 
Media (Wiley-Interscience, 1971); J. Jackson, Classical 
Electrodynamics (Wiley, New York, 1975). 
[11] Yu. Andreev et al, Phys. Lett. B 84, 247 (1979); 

D. McKay and J. Ralston Phys. Lett. B 167, 103 (1986); 
G. Frichter et al, Phys. Rev. Lett. 74, 1508 (1995), 
Erratum-ibid 77, 4107 (1996); R. Gandhi et al, Phys. 
Rev. D 58, 093009 (1998). 

[12] A. Donnachi and P. Landshoff, Phys. Lett. B 437, 
408 (1998); ibid B 518, 63 (2001); L. Bezrukov and 

E. Bugaev, Sov. J. Nucl. Phys. 33, 635 (1981). 

[13] S. Razzaque, J. Ralston, and P. Jain, in preparation. 
[14] S. Klein, Rev. Mod. Phys. 71, 1501 (1999), and references 
therein. 

[15] A. Cillis et al, Phys. Rev. D 59, 113012 (1999). 



